J V: ^ ' 
sF( K ? : f > " 


7 Oata ,7744'- X-621-72-234 

.: ! %7 v <; ' ?-&-r c yy 4y a a o- ? 44 oy j y 4 • : a 


" \ a/* 'S' 

^ Y - 


a v-4- u -'tv ■ _0 -''I; . *-'** ^ ii> '-' •’' "' 

""fy<l' y'- > J • Ait?". :■ 477.— _4 




IC -->a 


rv*- i , V > V- / S* Y. X •/■—<.. 1 -w r^^irtlsaS T V ' 7 f ' ••,'*• j p -- 

' ^ y * * •, %m . x ~"' \^ ■ - — - r --■ ->j y* - y, ^j y yi'. y s, . ' . -\ -y -* ■ \ 

A NUMERICAL STUDY ^§ 1^74 

iiA&Nsi©iM spliia^tAip^ 

p'-: , M^ : :J.r :.-. v 7*o -yYcyYyo 

aya ■ 4 - 4 , YV4 * - • Ty- a v - a- - ya yo ■ a-: ■- 


/ t j . ■>, 

! »v — > ■£« 7 . , a 


V- 4 . 


l^Yy 


>4, 


te>T . Y r ‘T 

a v -.’ 4 >: -a, 

.*MVY££ : 

’ JT 4 


Jy- -my^yuw 

7 ' ; < , a . ’ 4 a -s--' 

:> x..; v yy 


7 a - 4 Vyy-- v.a v y.o -, y : 

T-. 1 ^ - ■ \,-r-’ a ._'/>, . .y'S- 

• a . . v, t ^~s y ; ■ - 


t INO 1 1 I ' M IN Lf : VY I IN LT mUUCL 

7 -<• - >7 fe-' .: ■--. : 

Y4 >• - •_>"■ v />7, , a ■ W , ■ V’Sr^X'' 7 * 

r> x. - t. • •. • • > V/ l . vV- r- f ' O' • i\- N; 


>V^45V'^ 


’A' c- . 

— v> -■ yyy* \. 

H-:\2 


bi a Vw/-. ,.a r?-\ a a i r t ; a a r ^oa ■ v \ > 

^ a. “' a ' r S' ’■/ ■ - a- V s - " 

H .,<■ a.- v . • i * , ,-. ,y i ‘',4 a.//' s - - x . a : -x -r,.' s 

r- -oaV. V V J 'V v -V;^ ■-< : 5 ^ •• ; ^ r '^::o 

t'- :. .- ^ ^ ^-';a ;,: r ; • r-.l y > -r* ' 

;VaX. « - ■ 'Y^~ ■>•■'■; - o / ; > - ^ , vta^:,', *. ~o .O u ; '/ '-V 


'"• ■■ • fa ■■ -o ; 'o ,'«a _:, 75 f. ■< . O' •- f) o - jf- -'.v ay y;. a- , -:> v A'--v . 

v/r o - ^ vaa .’'VSr- :;>“.;v ; ,;? :-Y /-> a a^'- , - -a'/i Y .-• '/ V 

a: a r> ,04 ao, a y , § a a4 

. a> ;a«..o"'y-:axrva-?H.^oi^AHD>--'aafs c-vf^ 

J / .U' V, >yv , ^ ■ J - .. .„...'• •’••'.••O •• •_ - a.' -....' < — .• -a -A •'> '• ^ aa ; i , N a itz -:■•■■ y. r- ,-r 5 
a » > - <.< *.:«. • *. .<- ' ia a .. > ... ^ . _^- < ». J ,. — * , K ■ x, f-t- / -. 


‘#',aaa.aaa / 4o ■ ^ 

;••'•' '7 771 : 77 7 f| • ; TT: ,k : : r ■ : ’ vl > 5 ^ r 

!,v/ ■>fy i ' a? V - ".■<■'■ i ^ _ -, i - ■ N ” _■ ,u. ■~d=s} -a,- N . -■<• , -■( ■ *>v (... .y ■> O' x a.c • ’.y'.'..._r , a — , »j 

11 - ^a;r a'-lr^ 


•).'.• s) j."'.'"'^ a.'' y v a. •.',' ,■■^,.,,.'•, 0 ' a 

rs.V7' .' Xt 7'--a ^ v- >7--a T '- 

.•'? '^‘iLao.' -j ... i 1 ' - . -V ’ N t ,» N ' Vy 44 .' 4 : 

V a '. aai? ;4.\ a _ - •> .a-. ^ > '.: Vy'^a, 

o .' v ' \ 44 x ,-''' — ,"~ v>< . ' v !a ; ^"^yr- .xy :; f • V 


-YY.' - ^ vr YYcaY-.- 7^1 ■&\ a V ••> S-p? ; <Y w. 
aa. • : 4 /V Ya;' - ; ■■>>■ v ■ 77o 7a • a . T " ' a >o c '—a ■&$}/ • 7^- . . ,o 

a ; ■■-> -..a..'- u>o- '.'..'.y-aiaa i L'.¥ i 

ya :?...y aa ' -aa; a also- x.m.>- :<*&i 

t v\ „h i a'..-,... v^;--o T - ' s r;a .^a...l.-.' ^L’ik. o-r; . =Oj . .-.a,- .v-Tr 


11a: : 


c' . ^ * V ''|V \ 1"V -<S- >— ^ Y.‘ ' ” - '-yo' T ' '' ,Si Y* * ^ yl '7 . 7 } 7 aoY' ^vr« - ^ .- >-."»• -v , ' ■ , % v 

a-4 vf o-.-,.^. x june 1972 w^xYY-?a aa 4 

: .7 . 1 aa' Oil- ': ; ' : v! -.a- -■’ a-a. Ov ' ■-'•Vv 

' p^' *’ a 4 1 4/ •■■ j- VOo'-' - ’ si> 7 ' i-J w '/ ^ • ‘ V 4 --a7 O' 7 a^ 1 ; v ..'" a~ “■ (fS 

v v| .> , A -■ .\ ' >•- - - - Os.- ' a ' c - -yyy ■ --YY/ ,‘~ Y-. *•. ■ ri’ -'■••>• a- •■ ■ t r 


a - V. I 1 ■''“ /f"4 

» '■ : s' /..'■' y-' ■ 

ova 


/. 7 V I , ^ 7 77 ay ^; : \7y. 


.. <'. -'. lYaal - - V- w - ,x ■; 


' 11 a'^S-'V !77V' Y- -a74 0 ' Y:'la/' A 

>7Yw7/7V 1 ' r "■&*?..- ■ 77.-' / r --a ></’ ‘ ? a 5 7 <-.7' 

v r -.'a;a < it aa v- ..< -.-a :' ■•• -a . 'a-.-.. . - 


l-’Y. 


7; '" (a 

- ' < • f ; -a-y- 

y-y >•., a. 



-Tf.v ' ' /ro o ■■• 70 • . * \;-v 

7)4'„;-vy ; o- •- .;. 44 ; 

^lyy.! > yry-.j-pY^s^ . .yy 


■^17' yy^ 


' -fyj y 




y -7 >,v; 


GSFC 


:v y:'7a 


y o ' .->7 A ■ 


GODDARD SPACE FLIGHT CENTER 


*■ V V7/ T.. -, . 


; ~!\>.- > ,J 


GREENBELT, MARYLAND 


m. - a a 


7' Y 
K 


y 4 4:v4::ll>\’yl7l.a /; ; , '"t 
tvljl fl > 4 H-'l 77 - ■ (Yyl 1 4 >.a-y: S y • ^ylyyyylr^y .7 

> ^l< 4 4/4aa-';7 7;y 7/777. al a — > « , 


ffs. y v ’t>4 . 

Y, " 7 .V' • ya ; 


A /, 7L. ^y y • \ 1 y ^ s' n •• ■ y . 07 • - ! *v "O ^ - / iO 

yy. ;av a :,y... 7 .4 

,yy 4 b :x-6596 0 , s * B »s?s c sSp«Sic 
x!«!?“S: S SS“. S oSS , “ Voiland.y T al 


.-"a 1 ! ^ > 47 ^' 


N 72-28358 


'■ 7-7 
'/ ax. .-. 


IS” E : Eo : y“'™y^,s|f 


yi# 


4>:7xa'" 


Yl 71#: 


. or :i 'y 


.-O' ',7y.y.';a 

Ivl' Ja V 


;.c - . > r \ 

V 


X-621-72-234 


A NUMERICAL STUDY OF A THREEDIMENSIONAL SPHERICAL 
THERMOSPHERIC DENSITY AND WIND MODEL 


by 


H. Volland 

Astronomical Institutes 
University of Bonn 
Bonn, Germany 

and 


H. G. Mayr 

Goddard Space Flight Center 
Greenbelt, Md. USA 


June 1972 


Goddard Space Flight Center 
Greenbelt, Maryland 



PRECEDING PAGE BLANK NOT FILMED 


CONTENTS 

Page 

ABSTRACT vii 

1. INTRODUCTION 1 

2. BASIC EQUATIONS 2 

3. DEVELOPMENT IN TERMS OF TIDAL MODES 6 

4. THE VERTICAL STRUCTURE FUNCTION 8 

5. THE DETERMINATION OF THE EIGENFUNCTIONS 11 

6. MODE COUPLING 16 

7. THE FUNDAMENTAL SYMMETRIC DIURNAL TIDAL MODE (1,1,1) 18 

8. THE FUNDAMENTAL SYMMETRIC SEMI-DIURNAL TIDAL MODE 

(2,2,2) 30 

CONCLUSION 35 

REFERENCES 39 

ILLUSTRATIONS 

Figure Page 

1 Height profiles of mean temperature T Q , mean velocity of 

sound c 0 , mean collision number 0 and effective collision 
number v used in the model calculations 18 

e f f 

2 Solar EUV heat input Q i versus altitude used in the model 

calculations 20 

3 Magnitude (Fig. 3a) and phase (Fig. 3b) of the density coefficient 
p\ // 0 Q of the diurnal (1,1,1)- mode versus altitude generated 
by solar EUV heat input Q In ("EUV"), by the electric field of 
the Sq current ("EF") and by a heat source within the lower 


iii 


ILLUSTRATIONS— (continued) 


Figure 


3 

(cont) 


4 


5 


6 

7 


8 


8a 


8b 

8c 


9 


Page 

atmosphere ("TW"). The dashed lines give the sum of the 
three contributions. The thin solid lines have been calculated 
adopting heat sources Q : and Q n from Fig. 2 21 

Magnitude (Fig. 4a) and phase (Fig. 4b) of the meridional wind 

coefficient u^ 1 /c 0 of the diurnal (l,l,l)-mode versus altitude 

generated by the solar EUV heat input Q ni ("EUV") by the 

electric field of the Sq current ("EF") and by a heat source 

within the lower atmosphere ("TW"). The dashed lines give 

the sum of the three contributions 24 


Magnitudes (Fig. 5a) and phases (Fig. 5b) of the coefficients of 
pressure p \ /p 0 , temperature T J /T q , zonal wind v { a /c 0 and 
vertical wind w \ 1 /c Q of the diurnal (-l,l,l)-mode versus alti- 
tude generated by the sum of the three wave generators con- 
sidered 25 


Height functions determining the latitude structures of pressure 


and zonal wind of the (l,l,l)-mode 27 

Total pressure, meridional wind and zonal wind versus colati- 
tude d of the diurnal (l,l,l)-mode calculated for four different 
heights. Solid lines: magnitudes; dashed lines: phases 28 

Horizontal wind field of the (l,l,l)-mode as function of local 

time on the northern hemisphere 29 


Wind field at 100 km corresponding to the wind field of the 
Hough-function (1, -1) (Kato, 1966) . 29 

Wind field at 300 km due to solar EUV heat input 29 

Wind field at 300 km due to the sum of the EUV heat input and 

the contribution from the electric field of the S current. . . 29 

q 


Magnitude (Fig. 9a) and phase (Fig. 9b) of the density coefficient 
P 2 2>2 /p 0 of the semidiurnal (2,2,2)-mode versus altitude gener- 
ated by the solar EUV heat input ("EUV), by drag induced coup- 
ling ("DR") and by a heat source within the lower atmosphere 
("TW"). The dashed lines give the sum of the three contributions 33 


iv 



ILLUSTRATIONS- (continued) 


Figure 

10 


11 

12 


Magnitude (Fig. 10a) and phase (Fig. 10b) of the temperature 
coefficient T 2 2 /T q of the semidiurnal (2, 2, 2)- mode versus 
altitude generated by the solar EUV heat input ("EUV")> by 
drag induced coupling ("DR") and by a heat source within the 
lower atmosphere ("TW"). The dashed lines give the sum of 
the three contributions 

Magnitudes (Fig. 11a) and phases (Fig. lib) of the coefficients 
of pressure p 2 2 /p 0 , meridional wind u 2 2 /c Q , zonal wind 
v 2 2 /c 0 and vertical wind w 2 2 2 /c q of the semidiurnal (2,2,2)- 
mode versus altitude generated by the sum of the three wave 
generators considered 

Total pressure, meridional wind and zonal wind versus co- 
latitude 6 of the semidiurnal (2,2,2)-mode calculated for four 
different heights. Solid lines: magnitudes; dashed lines: 
phases 


Page 


34 


36 


37 


v 



PRECEDING PAGE BLANK NOT FILMED 


A NUMERICAL STUDY OF A THREE DIMENSIONAL SPHERICAL 
THERMOSPHERIC DENSITY AND WIND MODEL 

by 

H. Volland 

Astronomical Institutes 
University of Bonn 
Bonn, Germany 

and 

H. G. Mayr 
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Greenbelt, Maryland 


ABSTRACT 

Numerical calculations of the generation and propagation of the two impor- 
tant fundamental symmetric tidal wave modes - the diurnal mode (1, 1, 1) and 
the semidiurnal mode (2,2,2) — have been performed applying a realistic model 
thermosphere and taking into account heat conduction and the temporally and 
spatially varying ion-neutral collision number. Both wave modes are predomi- 
nantly generated by the solar EUV heat input. It is shown that the latitude struc- 
ture ofthe(l, 1, l)-mode which is identical with the Hough function (1, -1) within the 
lower non-dissipative atmosphere degenerates into the spherical function P at 
thermospheric heights. The pressure field of this mode constitutes the observed 
pressure bulge of the thermosphere the diurnal component of which peaks at 
15 h L. T. in the model of Jacchia and Slowey (1967). Its horizontal wind field 
agrees with the wind system derived by Kohl and King (1967) and by Geisler 
(1967). The electric polarization field of the geomagnetic Sq current generates 
a significant fraction of this wave mode at F layer heights. This wave compon- 
ent shifts the total horizontal wind system to earlier times by about 1 hour in 
agreement with ionospheric observations. The latitude structure of the (2,2,2)- 
mode is identical with the Hough function (2, 2) within the lower non-dissipative 
atmosphere. It degenerates to the spherical function P 2 2 at thermospheric 
heights. A component of this wave mode generated below 100 km can penetrate 
deep into the thermosphere and remains significant even at 400 km. Moreover 
an additional wave component is generated due to non-linear interaction between 
wave mode (1, 1, 1) and wave mode (2, 2, 2) via ion-neutral collisions at F layer 
heights. The interference between these two contributions and theXUV generated 
wave at thermospheric heights causes a phase delay between density and tem- 
perature of about 2 hours in the semidiurnal component which contributes to the 
observed delay between the times of the density and temperature maxima of the 
neutral gas at thermospheric altitudes. 
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A NUMERICAL STUDY OF A THREEDIMENSIONAL SPHERICAL 
THERMOSPHERIC DENSITY AND WIND MODEL 


INTRODUCTION 

The generation and propagation of tidal and planetary waves within the 
threedimensional thermosphere has been studied in a previous paper on an ap- 
proximate base (Volland and Mayr, 1971; referred to as paper I). The main 
assumptions made there were: 

a) Validity of perturbation theory 

b) Isothermal atmosphere 

c) Consideration of gravity waves only 

d) Constant dissipation factors. 

In this paper we want to extend this study by abandoning the points (b) to (d). 
We shall assume a thermosphere model with a realistic temperature profile. 

We shall introduce heat conduction waves besides the gravity waves and we shall 
take into account the temporal and spatial distribution of ion-neutral collisions 
in this model. We shall however again neglect the influence of viscosity waves 
and shall simulate their effect by an effective height dependent collision number. 
This study therefore will also be an extension to another previous paper which 
dealt with a two-dimensional thermospheric model in much the same way 
(Volland and Mayr, 1970 referred to as paper II). 

We shall do numerical calculations of the generation and propagation of 
two of the most import symmetric tidal waves at thermospheric heights — the 
fundamental diurnal (1,1, l)-mode and the fundamental semidiurnal (2,2,2)- 
mode. We shall study the influence of the solar EUV-heat input into the thermo- 
sphere above 100 km on the generation of the two tidal modes. Moreover, we 
shall estimate the signifiance of waves that are generated within the lower 
atmosphere below 100 km. 

Finally, we shall take into account waves generated by the electric polariza- 
tion field of the geomagnetic Sq current, by solar semi-diurnal gravitational 
forces and by non-linear interaction between the (1, 1, l)-mode and the (2, 2, 2)- 
mode resulting from the diurnal variation of the ion drag force. 
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The abbreviations used in (1) are the following: 

T 0 > P 0 , p Q , M Q mean values of temperature, pressure, density, and molecular 
mass 

T, p ,p , M perturbation amplitudes of temperature, pressure density 
and molecular mass 

v = (u, v, w) wind vector of the perturbation with its components north- 
wind u, westwind v, and vertical wind w 

— ♦ 

v hor = ( u > v ) horizontal component of v 

H 0 pressure scale height within the homogeneous isothermal slab 
v collision frequency of ion-neutral collisions 
u i » v i > w i components of the ion velocity 
9 ,k, r spherical coordinates 
s = sin 9 

c = cos 6 height above ground (r = earth’s radius) 
z = r - r 0 

fl angular frequency of the earth's rotation 
V coefficient of molecular viscosity 
K coefficient of molecular heat conductivity 
g gravitational acceleration 
c v specific heat at constant volumn 

Q heat input component which generates the perturbation 

$ gravitational tidal potential generating the perturbation. 

The collision number v is a function of time and space. Strictly speaking, 
it depends on the coupled dynamic system of neutrals and the ionospheric plasma. 
For convenience we assume v as an external known function of the form 


3 



( 2 ) 


V 


V 


0,0 0,0 


+ V. 


p +1^ p 

, 0 r 2, 0 r 1, 1 M, 1 


:os ( T - T ! !> 


where P n m (0) are the spherical functions in Neumann’s normalization, r is the 
local time, v nm (z) are height depending coefficients and r x x is the time of 
maximum to be specified below. 

The ion velocity will be determined from the Lorentz equation in which 
Coriolis force, inertia and pressure forces and the electron component are 
neglected: 


e - -* _ 
v in O'. - V) =_! (E + V X B) 
m. 


(3) 


This leads to the components (with w ~ 0) 


u i = - {(1 + S 2 s 2 ) u - 2 S c v + (1 + S 2 s 2 ) S u - 2 S 2 c v + 2 S 3 s c w } 


v i = — {2 §cu + v+ 2S 2 cu +Sv -S 2 s v; } 


(4) 


w i -i{2 8 2 scu + Ssv + 28 3 suc +S 2 sv„ + (l + 4S 2 c 2 ) 8 w } 

1 EE' ' E 


with 


A = 1 + 8 2 + 3 S 2 c 2 


s = sin 9 
c - cos 9 



v. 

i n 


1.83 x 10~ 12 — 





i n g/cm 3 


in g/mol 


4 



I 


n. 






B„ 


m. ion 


e i B o . 


m. 


ion gyro frequency- 


^e. ion 


mass 


charge 


v in ~ 1.56 p 0 /Mg /2 (Obayashi, 1964) 

ion-neutral collision number 

E = (E d , Ey E r ) electric field vector 


B 0 ~ 3 x 10 5 V s/m 2 magnetic induction of the earth’s co-axial 

dipole field at the equator. 

At heights above the dynamo region, the electric field is orthogonal to the geo- 
magnetic lines of force (w £ = - s u E /2c), and it is 8 > > 1. 

Thus, 

s 2 u - 2 c v_ 

U. ~ 

1 1 + 3 c 2 


v. ~ -5- at z > 150 km 

1 2 c 


2 s c u + s v 

£ 


Here, the wind induced component of the ion movement is primarily bound to the 
geomagnetic lines of force along which momentum can be transferred from the 
neutral to the ionized component. 
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Within the dynamo region it is S < < 1, and (4) becomes 


u. ~ u 
1 

v. ~ V at z < 130 km, (6) 


w A ~ 0 

and the ions essentially move with the neutral wind. 

A more sophisticated study has to take into account the electron component 
within the dynamo region (see e.g., Kato and Matsushita, 1969). However, since 
the ion-neutral drag becomes significant for the neutral gas dynamics m ainl y 
above the dynamo region we shall neglect those complications and s hall use 
Eq. (4) in the whole altitude range above 100 km. 


3. DEVELOPMENT IN TERMS OF TIDAL MODES 

It is well known from tidal theory that only a discrete number of eigenfunctions 
or tidal wave modes (the so-called Hough-functions) can exist within the atmos- 
phere (e.g., Chapman and Lindzen, 1970). These eigenfunctions will be written 
as (see Paper I) 


E n, m = ®n,m e X P (j ^ + j f fi t) (7) 

with 



n 


Here it is 


f = n ( “ +sJ V> 


a normalized frequency 


In this paper we introduced for convenience spherical functions in Neumann’s normalization. 
However, in paper I we used the spherical functions in Schmidt’s normalization. i It is 

. f i r „ £ 


@m, f 

n 


f ( n 


(n - m)! 


(n + m) ! 


and 


S m ' f , 

n , n 


: l/ ( n ’ + m )! 

1 (n’ - m) ! 


( n ~ m) 
(n + m) 


- S f , 

I n , n , m 
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w the angular frequency of the wave 
Ci the angular frequency of the earth’s rotation 
fi a the angular frequency of the earth's orbital rotation 
t universal time 

n, m, s integers specifying the type of the wave mode 
S 1 , coefficients of 0 to be determined below 

n , n , m n 

8 f = 1 . 

n , n , m 


The Hough-function of type (n, m, f) is determined from the eigenvalue 
problem of the two homogeneous equations (lb) and (lc) (setting E = $ = 0) 
together with 


div {v, } f 

hor J n,m 


j f p 


I '* 1 E f 


( 8 ) 


Here, h^ m is the so-called equivalent depth and p’ m is the pressure amplitude 
of the wave mode of type (n, m, f). 

The electric field E and the gravitational potential $ are considered as 
external forces, although strictly speaking, E in the case of the electric 
potential of the Sq current is generated itself by the neutral wind at E -layer 
heights. However, the energy necessary to generate E is only a small fraction 
of the kinetic energy of the driving wind which partly even originates below 100 
km altitude. Therefore, at E -layer heights this small fraction of the kinetic 
energy does not significantly influence the wind itself, and e can be calculated 
independently from the dynamo theory. However, E affects the neutral wind at 
F layer heights where the electromagnetic energy is here of the same order of 
magnitude as the kinetic energy of the wind. E, therefore, acts like a coupling 
link between E- and F -layer regions coupling a fraction of the kinetic energy 
of the wind at E-layer heights into the F-layer region. E tries to map the wind 
structure at E -layer heights into the upper part of the thermosphere. 
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We shall determine the eigenvalue problem in the next sections. Here we 
shall assume that the equivalent depth h f and the coefficients S f , of 

r n , m n , n , m 

©* are already known. Then, the perturbation amplitudes P, p, T, w and M as 
well as the heat input Q can be developed in terms of the wave modes E ^ : 


w 


M 


\ 

V 1 

r n , m 


Pn.m 

>=z: < 

jf 
n , m 

w f [ 

n , m , f 


M n m 
n , m 

J 

Q 

x n , m 

V- J 


E f 


( 9 ) 


Likewise the electric field components and the gravitational potential may be 
developed in terms of the eigenfunctions. We find then from (lb) and (lc) 


div v, 

• ho 


r " ‘ j n Z] | f \Z'- m + F + F g n.l 

n , m , f Pq n , m 1 


( 10 ) 


Here, F e f and F g f are functions of the components of the electric field 
and of the gravitational potential, respectively. They will be determined in the 
following sections. 


4. THE VERTICAL STRUCTURE FUNCTION 

Introducing (9) and (10) into (la), (Id), (le) and (If) allows to separate the 
various wave modes, and we obtain after elimination of p from (If) a system of 
four inhomogeneous linear ordinary differential equations of first order. They 
describe the altitude dependence of vertical wind w „ pressure p, f lim , tempera- 
ture T f and vertical heat flux F f = K d T f /dz of the perturbation ampli- 

n , m n,m n,m 

tudes of wave mode (n, m, f) which we write in concise matrix form 


j k n d z 


d g __ , 

— = Ke+h 


( 11 ) 
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with 




I \N f /C 

n , m 0 


Pn.n/Po 


T f /T 

1 n.m' 0 


\ F n >m /<^P0 C 0> 


(11a) 


K 


- j A 
0 
0 

\- 2 j A 


f (l-H 0 /O f 

j A - 2 j A 

0 j A 


- f 


f /k 


0 
0 

- 2 j k G 

- j A 


(lib) 


!i = / p7 r o 


2 ) »F-I„ 

0 

j T f / k - f F m f 

J J n,m' n , m 


f \ 

ra 


(lie) 


K Q‘ 

T f = n ' m 

n ' m n Po 


G = 2JL; v 0 

2 V n 0 


K g . 

"" — — f 

c p. 

p 


7-1 

7 


Fmf = !L- ; A = ^ ® - — - 
n,m 


2c 0 2k 0 H 0 


. Q 


y - ratio between the specific heats at constant pressure and constant 
volume 

c = ./ g 7 H q = velocity of sound. 
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The solution of (11) is straight forward and is fully described in paper II. 
Essentially, the procedure is to divide the real atmosphere into a number of 
homogeneous and isothermal slabs in which the elements of the matrices K and 
h are constant. Then one can transfer the physical parameters e in (11) into the 
characteristic waves: 


e = Q c 

where Q is a transformation matrix and 



( 12 ) 


(13) 


contains upgoing (a) and downgoing (b) waves of gravity wave type (i = 1 ) and heat 
conduction wave type (i = 2 ) the height dependence of which is 



exp (± j k Q q. z) 


(14) 


with the eigenvalues 


= f-A 2 - j f G + /- f 2 G 2 + 8 j f k G A 2 (15) 

The upper (negative) sign in the square root of (15) belongs to the gravity waves 
while the lower (positive) sign belongs to the heat conduction waves. This eigen- 
value has the same form as the eigenvalue of plane inertial waves (with A 2 , G » 
y , S 2 ) if we replace 


H 


/h f 


m 


= k 2 /k 2 7 f 2 = S 2 /y 


(16) 
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where k x is the horizontal wave number of the plane waves (Pitteway and Hines, 
1963). However, contrary to plane waves, the tidal waves have only a discrete 
number of (generally complex) horizontal wave numbers. 

Such transformation from the physical parameters into the characteristic 
waves as well as the slab method have the advantage that analytic solutions exist 
appropriate for numerical calculations and that the radiation condition can be 
adopted as boundary condition (see paper II). 

The inhomogeneous term h in (11) is responsible for the excitation of the 
wave mode. It contains the external forces due to heat input J f , electric fields 
Fe and tidal gravitational potential F g n m . Since our system (11) is linear, 
particular solutions of (11) can be superimposed to give the complete solution. 
Thus, we are able to separate the influence of the various excitation mechanisms 
by setting all but one of the elements J, Fe and Fg equal to zero. 

There remains the additional term Fm f in (11c) which describes the 
temporal and spatial variation of the mean molecular mass. It is only a fictive 
external source because this variation is generated by the dynamic processes 
within the thermosphere and therefore is originally generated by the real external 
forces. As it is well known, the thermosphere is neither completely mixed nor 
in a complete diffusion equilibrium. An exact calculation of the thermosphere 
dynamics including diffusion processes must therefore be based upon the coupled 
system of hydrodynamic and diffusion equations for the various constituents. 

In the following calculations we shall neglect the term F m f because it 
can be shown that F m does not influence significantly the pressure and density 
amplitude of the tidal waves. 

It follows then from (If) that the magnitude of the temperature becomes in- 
creased by Fm. Therefore, in the following results our relative temperature 
amplitudes may be underestimated as the result of the neglect of the variation 
of the molecular mass. 


5. THE DETERMINATION OF THE EIGENFUNCTIONS 

We now return to the remaining problem of determining the eigenfunctions 
and the equivalent depths of the tidal modes. In paper I, we used an approximate 
way to calculate the coefficients S l n _ n > _ m of the eigenfunctions 9 f n<m and the 
equivalent depths h nm . The assumptions made were the following: 
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a) Replacement of the viscosity force by an effective Rayleigh friction 
parameter: 



where Z v . s is a positive height denpendent number. 

b) Collision number v only height dependent 

c) Zero ion velocity. 

Here, point (a) is the most restrictive one because it is equivalent to an 
outfiltering of two kinds of characteristic waves, the so-called viscosity waves, 
which mainly influence the horizontal velocity. We shall retain point (a) because 
otherwise the mathematical difficulties would become overwelming. Moreover, 
it can be shown that the number Z vis is small when compared with the corre- 
sponding ion drag term 


Zcoi = */n 


(18) 


below 300 km. Therefore, viscosity has in fact not much influence on thermo- 
sphere dynamics in that height range (Geisler, 1967; Kohl and King, 1967). At 
greater heights, however, the viscosity tends to keep the horizontal velocity 
constant with altitude. It is possible to simulate this influence by an effective 
collision number above 300 km altitude (see paper II). 

We want now to abandon points (b) and (c) in our treatment. We shall show 
explicitly the way to do this for the fundamental symmetric tidal mode (1, 1, 1) 
which is the trapped mode (1, -1) within the lower atmosphere according to 
Kato (1966). The treatment of the other tidal modes follows the same line. 

In the case that points (a) to (c) were applied we found in paper I the solutions 
for pressure amplitude p and horizontal winds u and v of wave mode (1, 1, 1) 

(n = 1; m = 1; f = 1) from Eqs. (lb), (lc) and (17) 


= <Pl.I P l.l+P5.1 P 3.1> e,T 


u = u! 


1,1 sin 0 


(19) 
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v = (v 


1 

1, 1 


1, 1 


+ V 


p ) 

3 , 1 3 , 1 ' 


eJ T 
s i n 8 


with 



( f k + 2/3) j x x ( f k - 14/15) _Cg_ i 

3y£A p 0 Pl>1 ’ 1,1 y^A p 0 Pl 


P 


l 

3 . 1 


8 i 

45 A Pl - r 


3 , 1 


4 

15 (f k - 14/15) 


v 


l . 
l, l’ 


A= f k - 4 V 3 - 4 / 15 


€ = — r; f k = 1 - j (Z vis + Z col ); T = at+^ (local time). 
c o 

From Eq. (10) it follows according to paper I (Eq. Il/7) approximately 



a 2 r 2 A 

2 g (f k - 1/3)’ 


( 20 ) 


We retain in the following the general form of (19), however with coefficients 
to be determined below. Moreover, we introduce 


V (u - u.) 

V (V - V.) 


n 


y y * (A| t 1 u l, 1 + ®r,l V l, 1 + 3 r , 1^ ^r, 


sin 


« e 


Jfr ( 21 ) 


f = 0 


with 


q = 


l r s l 

for 

1 r = 0 • 
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Here, for the meridional component u, it is r = 2, 4, 6 . . . , and for the zonal 
component v, it is r = 1, 3, 5 . . . 

The coefficients A | > t and B \ j can be derived from (2) and (4) applying 
well known theorems of spherical functions. For example it is 


A 


i 

r , 


1 


(2 r + 1) 
r (r + 1) n 



0, 0 


- V 2 , 0 /2 + 3 v 2,0 c2/2 ) 

(1 + S 2 + 3 §2 c2) 


P 


2. 


P r , 1 ^0 


"4 8 2 


.-28c 


for 


r = 2 n' 
r = 2 n' + 1 


( 22 ) 


Similar formulas exist for the coefficients B \ 1 . The integration of (22) is 
straightforward . 

The electric polarization fields in (4) give rise to the coefficients 


d 1 , - (2 r + 1) 

r , 1 “ 


r (r + 1) f) J 


e_iT ( V 0.0 - ^2.0 /2 + 3 v 2. 0 c2/ 2) s P r. x d C < 


2 C V E 

1 + 3 c 2 


U E 

2~c 


(2 n' 

: o r n = J 


and z > 150 km 


(23) 


2 n' + 1 


dj x ~ 0 for z < 150 km. 

The last assumption is reasonable because, below 150 km, it is dj j 00 8 Q . 
8 and v 0 0 are both small in this height range. 

In (23), the electric polarization field is due to the Sq current and has the 
form (Volland, 1971) 


14 



( 24 ) 


u £ = - 3 c (1 - 2 c 2 ) j V e e jT 
v E = - (1 + 2 c 2 ) V e e jT 


with 


V e = 40 m/sec during moderate solar activity 
c = cos d\ s = sin 9; c Q = velocity of sound 
r local time, zero at midnight. 

Since the symmetric diurnal component of the solar gravitational tidal po- 
tential is zero (Chapman and Lindzen, 1970), we set furthermore 

j = 0- (25) 

For the moment we shall neglect the series with the periods 0 and 2r on the 
right hand side in (21). We shall discuss their influence later. Introducing now 
(17), (19), (21), and (25) into (lb) and (1c), applying well known formulas of 
spherical harmonics and truncating the development into series of spherical 
harmonics after the terms with n = 4, we obtain from (lb) and (lc) four linear 
inhomogeneous equations for the coefficients Uj 1 t , v* t , v 3 t and p 3 t with the 
solutions 


3 , 1 

> 5.1 


1 = a. p j j + b. V e (i = 1, 2, 3, 4) 


(26) 


Eliminating the wind components of (26) in (10) and applying the same approxi- 
mation as in paper I (Eq. II/7), we obtain the equivalent depth and the external 
forces 


h 


l 

i . l 


2 ^Hq c o 

3 p Q (j a t + a 2 + a 3 ) 
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F S I , i ~ 0 


F e 1 
r l.l 


3 V 

(jb.+b.i b,)JL 

? 0 for z 

0 



(27) 


Evidently, the accuracy of this method depends on the degree of truncation 
of the series of spherical harmonics and can in principle be improved by taking 
into account higher order terms in the series. 


6. MODE COUPLING 

In our treatment, we use perturbation theory and neglect products of the 
wave parameters as small values. Thus, we eliminate any coupling between 
modes of different wave domain numbers (n, m, f). An estimate about the contri- 
butions of the non-linear terms led to magnitudes not greater than 20% of the 
linear terms (see paper I). Such errors should be tolerable in view of the com- 
plex problem. In spite of perturbation theory, we expect two other coupling 
mechanisms between the various modes. The first one is due to the diurnal varia- 
tion of the collision number v in (4). The second one is related to the change 
with height of the latitude structure of the eigenfunctions. 

We first consider the effect of collisions. In (21) we already added two series 
with periods zero and 2 r which followed from products of the kind 


p i,i *1.1 cos ( T - T U )Pi t i 


i,i' 


cos (r - t j ) 


(28) 


and equivalent expressions for the zonal wind component v j . The determina- 
tion of the coefficients in (21) was outlined in section 5. Evidently, the coefficients 
A° j , A J j etc. lead to mode coupling from the (1,1, l)-mode into the semi- 
diurnal (n, 2, 2)-modes and the planetary (n, 0, 0)-modes with periods 2 r and 
zero, respectively. This coupling does not affect the (1, 1, l)-mode itself. 


On the other hand, the modes (n, 2, 2) and (n, 0, 0) couple part of their energy 
into the (1, 1, l)-mode via neutral collisions. From the observations of thermo- 
spheric density variations, it is obvious (Jacchia and Slowey, 1967) that in the 
case of the symmetric modes only the modes (2, 0, 0) and (2, 2, 2) have significant 
amplitudes besides the dominant (1,1, l)-mode. The degree of coupling can be 
estimated from the ratios 
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"l.lK.O 


V l.l (^ 1.2 


( 29 ) 


and 

2 v o,o ^ v o,o ( u )i,i 

where (u) f are the maximum amplitudes of the meridional winds of the respec- 
tive modes , m and the factors 2 stand for the splitting of products of the kind (28) 
into modes of period r and 3 t . 

Assuming v ~ v l , we estimate from the result of paper I ratios of 
(29) of 0.05 and 0ll5, respectively. They are smaller than the contributions 
from the non-linear terms of the wave parameters. Thus, it is justified to 
neglect this coupling into the (1, 1, l)-mode due to ion-neutral collisions. On the 
other hand, the coupling from the (1,1, l)-mode into the modes (2, 0, 0) and (2, 2, 2) 
expressed in the two series in (21) is by no means negligibly small. It acts like 
an "external" force for the generation of these two modes. 

The second kind of mode coupling occurs in connection with the height change 
of the latitude structure of the eigenfunctions as we see this for example from 
the coefficient 


a ° (30) 

4 45 (fjj - 4 f k /3 - 4/15) 

of the (1, 1, l)-mode (see Eq. (19) and (26)). Even in the simplest case of (19) 
is f k a function of altitude. Within the lower atmosphere, it is f k =1 and a 4 = 
0.3. Above about 200 km, it is | f k | > > 1, and | a 4 | < 0.05. It is the altitude 

range between 100 and 200 km, where the change of the latitude structure of the 
modes is most drastical (see Fig. 6). 

This altitude dependence of the mode structure is however in contradiction 
to the method of the separation of the variables applied in this paper which pre- 
supposes the eigenfunctions m to be altitude independent. The effect of this 
altitude dependence is that, at "any height, part of the energy of the modes is 
either partially refelcted or transferred into an other mode of the same wave 
domain number n. Unfortunately, this mode coupling is a completely open 
problem and not yet solved. One can perhaps argue that because of the observed 
insignificant contributions of the modes with wave domain numbers n > 2 at 
thermospheric heights the degree of coupling from the higher order modes into 
modes of lowest wave domain numbers should be very small, whereas on the 
other hand the thermosphere filters out the higher order modes which are 
generated due to this coupling from the lowest order modes. We shall not try 
to estimate the degree of coupling in this paper but rather leave it to further 
studies. 
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7. THE FUNDAMENTAL SYMMETRIC DIURNAL TIDAL MODE (1,1,1) 


We are now prepared to do numerical calculations for the generation and 
the propagation of tidal waves at thermospheric altitudes. We shall restrict 
ourselves in this paper to the two important symmetric modes: the fundamental 
diurnal tidal (l,l,l)-mode and the fundamental semi-diurnal mode (2,2,2)-mode. 
In this section we start with the diurnal mode (1,1,1). 

We select as in paper I a height profile of the mean collision number 
plotted in Figure 1 as solid line which is roughly proportional to the electron 


c Q (km/ s) 

.2 1 



Figure 1. Height profiles of mean temperature Tq, mean velocity of sound Cg, mean collision 
number Vq q and effective collision number v e ff used in the model calculations. 
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density. Above 300 km we add to this collision number a Rayleigh friction term 
v = Qz (see Eq. (17)) which simulates the influence of molecular viscosity 
The sum v e = v + v is constant above 300 km (the dashed line in Fig. 

e f f 0,0 vis 

i). 

The spatial and temporal distribution of the collision number is given by 
Eq. (2) with the coefficients 


u 2,0 ® 7 6 V o , 0 


v - V 

1,1 0,0 


Tj x = 12 h L.T. 


The mean temperature T Q , the mean pressure p q and the mean molecular 
mass M Q are taken from the Jacchia model (Jacchia, 1965) for an exospheric 
temperature of To = 1000°K. The temperature and the velocity of sound for 
this model are plotted in Figure 1 versus altitude. 

As solar heat input we adopted three height profiles: 


Q oc p 


Q IX a exp (- z/H); (H = 50 km) 


(31) 


z ^ 160 km 


Qni 


for 


exp (- z/H); (H = 27.5 km) z < 160 km 


V 


Below 100 km the heat input is taken as zero. In all three cases maximum 
heat input occurs at 12.00 h local time. The heat input Q was used in paper I, 
the heat input Q XI was used in paper II. 
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The height profiles of these three heat inputs are plotted in Figure 2. The 
numbers E j give the total heat input per volume and time above 100 km: 


Ej = Q. d z (erg/cm 2 sec). (32) 

z g = 100 km 

These numbers are chosen such that the corresponding density amplitudes at 
400 km altitude agree with the observed data. 

Figure 3 show the coefficient p\ l / P 0 of the relative density amplitude in 
magnitude (Fig. 3a) and phase (Fig. 3b) versus altitude generated by the respective 



-10 -9 -8 -7 -6 -5 -4 

log Q (erg/ cm 3 s ) 


Figure 2. Solar EUV heat input Q; versus altitude used in the model calculations. 
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Figure 3. Magnitude (Fig. 3a) and phase (Fig. 3b) of the density coefficient p\ ^ Pq of the 
diurnal (1, '1, l)-mode versus altitude generated by solar EUV heat input Qm (“EUV"), by 
the electric field of the Sq current (“EF”) and by a heat source within the lower atmosphere 
(“TW"). The dashed lines give the sum of the three contributions. The thin solid lines 
have been calculated adopting heat sources and from Fig. 2. 







three heat sources Q j, Q IX andQ in (Symbols: I, II and III). Also plotted in 
Figure 3 is the density of a free internal upward propagating wave with magni- 
tude IpJ J / Pq = 0.1 and time of maximum of 18.00 h L.T. at 100 km altitude 
(Symbol:’ T.W). Finally, the density amplitude due to the electric polarization 
fields of the geomagnetic Sq current generated by the term Fe j ^ in (27) is plotted 
in Fig. 3. (Symbol: "EF"). From Fig. 3a we notice that the wave from below 
decays rapidly between 100 and 200 km and becomes insignificant above about 
180 km altitude. This is contrary to the result in paper II where we over- 
estimated the influence of the wave from below due to the neglect of the Coriolis 
force. However, above 230 km the wave from below although of insignificant 
amplitude again increases slowly with height which is in agreement with the 
result of the two dimensional model in paper II. It shows the negligible influence 
of the Coriolis force on wave propagation in that height range. The density 
amplitude excited by the electric Sq field is only a small fraction of the total 
density amplitude in the entire altitude range considered and generally can be 
neglected. In the case of the wave generated by the solar XUV-heat input, we 
know from the geomagnetic Sq variations that the corresponding maximum hor- 
izontal wind between 100 and 200 km is of the order of 100 m/sec and the relative 
density amplitude of the order of 0.05. (Kato, 1956, 1966). The same order of 
magnitude of the density variations is determined from rocket born and satellite 
drag observations (Spencer et al., 1966; Taeusch et al., 1968; King-Hele and 
Hingston, 1967, 1968). These facts suggests that heat source III may give the 
best fit to the observations. The total density amplitude is the sum of the com- 
ponent generated by the XUV-heat input Q XII , the electric field induced compo- 
nent and the wave from below. It is plotted as thick dashed line in Fig. 3. 
Evidently, the wave from below predominates within the first 50 km above our 
lower boundary of the model. 

The total heat energy of heat source III above 100 km is E m = 3.3 erg/cm 
sec. This value is at the upper limit of the results of Hinteregger et al. (1965) 
from observations of the solar spectrum which is another reason to exclude 
heat source I. According to Hays (1970) the effective solar heat input is roughly 
proportional to the mean pressure. Our calculations show that this appears to 
be a reasonable approximation above E-layer heights. Within the E-region 
however it seems to be too high. Our normalized heat input J j 1 introduced in 
(11c) is 


jJ >x = 0.17 at z ^ 160 km (33) 

This value is greater by a factor of at least 40 when compared with the corre- 
sponding number within the Ozone layer at stratospheric heights (Chapman and 
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Lindzen, 1970). At 100 km, this value reduces by a factor of 10 to J t = 0.017 
which is still greater by a factor of 4 when compared with the Ozone layer 
source. 

The amplitude of the wave from below in Fig. 3a corresponds to a heat in- 
put of J j = 0.017 within the entire lower nondissipative atmosphere below 
100 km. Because the estimated energy input seems much smaller at strato- 
spheric heights, this may be an overestimate. However, since the penetration 
depth of this evanescent wave is limited to a few scale heights, only the heat 
input into the mesosphere and into the lower thermosphere accounts actually 
for the amplitude of this wave at 100 km altitude, and the exact number of this 
heat input is not known. 

In Fig. 3b we plotted the time of maximum of the density versus altitude. 

The XUV-generated wave peaks between 14 and 15 hours local time above 250 
km altitude which is in agreement with the satellite drag data (e.g. Jacchia 
and Slowey, 1967). Below 150 km where the wave from below predominates the 
time of maximum shift toward 18 hours local time. A density maximum between 
16 and 18 hour is just necessary to explain the phase of the Sq current (Tarpley, 
1970; Volland, 1971). The electric field induced density component peaks near 
22 hours local time above 180 km and shifts the total density maximum to 15 
hours at 400 km altitude. 

In Fig. 4a we plotted versus altitude the magnitudes of the normalized 
meridional wind components uj /c Q which are generated by the EUV heat in- 
put (heat source III in (31)), by the electric field of the Sq current and by the 
wave from below. The dashed line gives the sum of all three contributions. 

Here the electric field induced wave has the same order of magnitude as the 
EUV-induced wave at F2 layer heights. The wave from below dominates below 
130 km and becomes insignificant above 160 km. The reason for the difference 
between the large wind amplitude and the small amplitude of density (and there- 
fore pressure, temperature and vertical wind) of the electric field induced wave 
(see Fig. 3a) is that the electric field induced component of the horizontal wind 
field is nearly divergence free which implies a very small driving term Fe for 
the pressure in (27). 

The phases of the meridional winds shown in Fig. 4b vary between about 
21 h and 5 h local time above 120 km. The phase of the total wind remains be- 
tween 0 h and 2 h local time in the height range from 100 to 400 km. This is 
just the phase necessary to explain the phase of the geomagnetic Sq current 
(Tarpley, 1970). 

In Fig. 5 we plotted magnitudes and phases of the relative amplitudes of 
pressure p* 1 /v 0 , temperature T j j/T 0 , vertical wind w* 1 /c Q and zonal wind 
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temperature T] ^/Tq/ zonal wind vj j/cq and vertical wind wj ^/cq of the diurnal (1, 1, 1)« 
mode versus altitudes generated by the sum of the three wave generators considered. 





V 1 , i / c 0 generated by the sum of EUV heat source III, electric field and wave 
from below. This figure shows a general behavior expected from the simplified 
theory in paper I. Pressure, vertical wind and zonal wind tend to increase with 
height to an asymptotic value. The temperature remains nearly constant above 
200 km. 

An interesting feature in Fig. 5a is the relative minimum of temperature 
and pressure near 120 km which is due to the steep mean vertical temperature 
gradient in this height range. The exact shape of the horizontal wind profiles 
above 300 km in Fig. 4a and 5a is subject the some ambiguity in our model be- 
cause we introduced an effective collision number in order to simulate the in- 
fluence of molecular viscosity. However any realistic change of this effective 
collision number above 300 km although altering somewhat the horizontal wind 
profile does not influence very much the temperature and pressure profile in 
Fig. 5a. 

The times of maximum in Fig. 5b are nearly constant with altitude above 
200 km indicating the quasi-evanescent nature of this wave mode at thermospheric 
altitudes. 

Figure 6 gives the magnitudes of the coefficients a 3 /a 2 , and b /b from 
Eq. (26) which are responsible for the latitude structure of pressure (and there- 
fore temperature, density, and vertical wind) and zonal wind. The number a 
which is identical with the coefficient S x 3 x in (7a) decreases from 0.3 at 1$0 
km to a value of | a 4 | <0.02 above 250 km. Thus, the influence of the spherical 
harmonic P 3 , in (19) decreases with height, so that above about 250 km only 
the spherical harmonic P x t determines the latitude structure of the pressure. 

A similar situation exists for the zonal wind which is generated by EUV heat. 

Here the ratio I v 3> j/v/ j I decreases from a factor of four at 100 km to a value 
<0.2 above 200 km altitude again diminishing the influence of the spherical 
function P 3 x in (19), so that the zonal wind becomes nearly independent of 
latitude. However the significant amplitude of the zonal wind structure generated 
by the electric field is mainly determined by the ratio |b 3 /b 2 | which remains 
nearly constant above 200 km. 

In Fig. 7 we show the latitude structures of relative pressure p/p Q , meri- 
dional wind u and zonal wind v of tical mode (1,1,1) for four different heights. 

Here we notice again the degeneration of the pressure function from the Hough- 
function (1, -1) at 100 km to the spherical function P x x = sin 8 above 200 km 
altitude where it is identical with Jacchia's pressure bulge (Jacchia and Slowey, 
1967). The meridional wind does not change its structure in our theory because 
of the assumptions made in (17). At 100 km the zonal wind becomes zero at 25° 
latitude and changes its sign there. This results from the influence of the 
spherical function P 3 x . The abrupt phase transition becomes more gradual 
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Figure 6- Height functions determining the latitude structures of pressure and zonal wind 

of the (1, 1, l)-mode. 

with increasing height although it is still evident at 300 km altitude. Since for 
the EUV induced wave the phase transition is already small above 300 km (see 
Fig. 6 and paper I), it is the electric field induced wave which is responsible 
for the phase change in Fig. 7c. The electric field induced wave essentially 
maps the zonal wind structure of the E-region into the F2 region. 

The horizontal wind structure is plotted versus local time on the northern 
hemisphere in Fig. 8. In Fig. 8a the wind at 100 km is identical with the wind 
field of the tidal (1, -1) mode (Kato, 1966). The maximum northerly wind blows 
at 0 h local time at the poles. Moreover, we observe the reversal of the zonal 
wind direction at 25° latitude. In Fig. 8b we plotted the wind field of the EVU 
generated wave at 300 km altitude. Here the reversal of the zonal wind direction 
has disappeared, and the maximum northerly wind blows at 2.4 h at the poles. 

Fig. 8c finally gives the wind field at 300 km due to the sum of the EUV-heat, 
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Figure 7. Total pressure, meridional wind and zonal wind versus colatitude i9of the diurnal 
(1, 1, l)-mode calculated for four different heights. Solid lines: magnitudes; dashed lines: 
phases. 

electric field and wave from below. Now a reversal of the zonal wind direction 
again exists near 30° latitude in the morning and the evening hours, and the 
maximum northerly wind has shifted to 1.2 hours local time. Besides this phase 
shift, both wind fields in Figs. 8b and 8c have similar structures apart from 
the morning and evening sectors at low latitudes. 

Figure 8b is similar to the wind fields derived from the Jacchia model by 
Kohl and King (1967) and by Geisler (1967). (Note that we consider only the di- 
urnal component whereas, in Jacchia's model a small but significant semidiurnal 
component is included). Kohl et al., (1968) have calculated from the ionospheric 
continuity equation diurnal variations of the critical F2 frequency and of the 
critical height for various stations. In these calculations the horizontal neutral 
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Figure 8. Horizontal wind field of the (1, 1, 1)-mode as function of local time on the 

northern hemisphere. 



winds which cause vertical ionization drift have been adopted from the model of 
Kohl and King (1967). From a comparison between the experimental data and 
their theory, Kohl et al., (1968) concluded that a wind model with its phase ad- 
vanced by 1.5 hours would improve the data fit very much. They already sug- 
gested electric fields as a possible cause for this shift. Our result in Fig. 8c 
verifies this suggestion. The electric field of the Sq current indeed prevents 
the phase of the F2 layer winds to shift by more than 1 .2 hours with respect to 
the wind at E-layer heights. 


8. THE FUNDAMENTAL SYMMETRIC SEMIDIURNAL TIDAL MODE (2,2,2) 

In order to calculate the generation and propagation of the symmetric semi- 
diurnal mode (2,2,2) we used the same model of the mean thermospheric proper- 
ties and of ion drag as for the mode (1,1,1) in section 7. We furthermore adopted 
heat source Q n i in (31). However we reduced the total heat input of this source 
by a factor of 5/32 to E m = 0.52 erg/cm 2 sec. The factor 5/32 corresponds 
roughly to the ratio between the semidiurnal and the diurnal heat input compo- 
nents in a development into spherical functions (see paper}, Eq. (1/9)). Because 
of the apparaently small semidiurnal component of the tidal winds which generate 
the geomagnetic Sq current (Kato, 1956; Stening, 1969; Tarpley, 1970) we expect 
small effective electric polarization fields of this period at F2 layer heights. 
However, a solar semidiurnal gravitational potential of 


0 = - 3.7 x 10 -7 P 2 2 cos 2 r (km 2 /sec 2 ) 


(34) 


(Siebert, 1961) acts like an additional driving mechanism for such a wave. In 
the case of the nondissipative lower atmosphere it can be shown (see paper I 
Eq. (H/34)) that the gravitational driving force in (11c) becomes 


l(F g)2, 2 l ^ 


35 

6 U 2 r 2 



- 10 -5 


(35) 


where 4> 2 2 is the coefficient of $ in a development into a series of the Hough- 
functions. The number in (35) must be compared with the corresponding heat 
generator in (11c) which is 


30 



T2 

J 2, 


0.005 x 
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at stratospheric altitudes. At thermospheric heights it is according to paper I, 
Eq. (11/26) 


(F g) 2 J 


6 * 1 . 2 

f k fi 2 r* 


~ 2 x 10' 6 


while J 2 - 3 x 10~ 3 there (see (33)). Therefore the tidal gravitational force 
does not significantly contribute to the semidiurnal solar tides at thermospheric 
heights. 

However, the non-linear coupling from the tidal (l,l,l)-mode into the (2,2,2)- 
component via the diurnal term l of the collision number (see Eqs. (2) and 
(21)) may give rise to significant amplitudes. The coefficients A r> j and B r l in 
(21) have been calculated according to the scheme given in (22), and the wind 
coefficients u 1 l and v} t from Figs. 4 and 5 have been adopted in order to de- 
termine a virtual external force in (10) or (11c) which acts like the generator of 
semidiurnal waves in addition to the EUV- solar heat input. 

Finally we added to these two wave generators a wave generator situated 
within the lower atmosphere below 100 km. According to the few direct obser- 
vations at E-layer heights (Greenhow and Neufeld, 1961; Rosenberg and Edwards, 
1964; Roper and Elford, 1965; Zimmer mann and Marcos, 1967; Murphy and Bull, 
1968; Bernard, 1971) the tidal (2,4,2) component seems to dominate at E-layer 
heights with wind amplitudes of the order of lOm/sec. This wave has disap- 
peared above 200 km altitude which can readily be explained by the filtering ef- 
fect of the thermosphere dynamic system (see paper I, Eq. (1/ 49)). At higher 
altitudes the (2,2,2) component is the dominant semidiurnal component (see 
paper I, Eq. (1/14)). Due to our ignorance about the semidiurnal (2,2,2) com- 
ponent at E-layer heights we can do only sample calculations of the wave from 
below such that the numerical results are consistent with the observations 
above 200 km. For this reason we selected a density amplitude at 100 km for 
the wave from below of p 2 2 /p = 0.008 peaking at 2.7 hours local time. These 
two numbers are not completely arbitrary because they correspond to a heat 
input of J 2 = 1.4 x 10“ 3 with maximum heat at 0 hour local time within the 
height region between 70 and 100 km where this wave is a trapped wave. This 
is equivalent to a solar heat input proportional to cos x (x is the zenith distance) 
into this height range with maximum amplitude of J = 0.004 at 12 hour local time. 
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Figure 9a shows the relative density amplitude versus height generated by 
the three wave generators considered as well as the sum of the three contribu- 
tions (dashed line). The (2,2,2)-mode being a trapped mode between 70 and 100 
km becomes a propagation mode above 100 km altitude. Therefore between 100 
and 120 km we notice in Fig. 9a a nearly exponential increase of the tidal wave 
from below (Symbol "TW"). Above 120 km however, this increase is compen- 
sated by the dissipative effects of heat conduction and ion drag transforming the 
(2,2,2)-mode into a quasi-trapped mode with nearly constant amplitude with 
height. The behavior of the EUV-generated wave (Symbol "EUV") is quite normal 
and similar to that of the diurnal mode in Fig. 3a with the difference that its 
amplitude is already more pronounced below 200 km due to the propagation 
characteristics of this semidiurnal mode within the lower thermosphere. The 
drag induced component (Symbol "DR") which resulted from the non-linear 
coupling between the diurnal and semidiurnal mode is unimportant below 250 km. 
However, it becomes a significant fraction of the total wave above this height. 

The wave from below dominates below 200 km and is still non-negligible above 
this height. 

Figure 10a gives the relative temperature amplitudes versus height of the 
three components with a similar behavior as in Fig. 9a apart from the fact that 
the wave from below ceases to dominate already above 160 km. The maximum 
times of density and temperature are given in Figs. 9b and 10b. The EUV- 
generated wave peaks near lh L.T. above 200 km altitude while the drag induced 
components peak near 3h L.T. The density amplitude of the wave from below 
changes its phase with height from 2.7h at 100 km to 7.4h at 400 km. The tem- 
perature amplitude of the wave from below is phase advanced by about 3 hours 
with respect to the density. This phase advance causes a time delay of 1 to 4 
hours between temperature and density of the total wave (the dashed lines) be- 
tween 200 and 300 km. This time delay in the semidiurnal component is con- 
sistent with the results from a Fourier analysis of the observed neutral density 
and ion temperature (Volland, 1970) and contributes to explain the difference be- 
tween the times for the density and temperature maxima. 

At 400 km height the ratio between semidiurnal and diurnal density is ac- 
cording to Figs. 3a and 9a 


I P 2 _ 2 />° i . 1 1 ~ 0-08- 

This number agrees well with the corresponding number from the Jacchia-model 
(Jacchia and Slowey, 1967) which is (see paper I, Eq. (1/14)) 
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Figure 10. Magnitude (Fig. 10a) and phase (Fig. 10b) of the temperature coefficient Tj 2^0 
of the semidiurnal (2, 2, 2)-mode versus altitude generated by the solar EUV heat input 
(“EUV"), by drag induced coupling (“DR”) and by a heat source within the lower atmosphere 
("TW”). The dashed lines give the sum of the three contributions. 
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0 . 12/12 
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(Note that Jacchla’s exospheric temperature is merely a measure of the ob- 
served density). 

Figure 11 shows the coefficients of horizontal winds, vertical wind and 
pressure generated by the sum of the three wave generators considered. Since 
the drag induced fraction of the horizontal winds exceeds the EUV-induced winds 
between 200 and 300 km and since the phases of both components differ by about 6 
hours the amplitudes and phases of the total horizontal winds are rather variable 
with height. Above 300 km this may not be realiatic because of our simplifying 
assumption about viscosity. However pressure and vertical wind behave in a normal 
manner. An important feature in Fig. lib is the strong phase change with height of 
the horizontal winds which covers about 2/3 of a period between 100 and 300 km. 

It is this phase transition which prevents the semidiurnal mode from gen- 
erating a significant contribution to the geomagnetic Sq current (Tarpley, 1970). 

Figure 12 finally shows the change with altitude of the latitude structure of 
relative pressure, meridional wind and zonal wind of the semidiurnal (2,2,2)- 
mode. Here again we notice the transition of the pressure structure from the 
Hough-function (2,2) at 100 km to the spherical function P 2 2 above 200 km (see 
paper I). Moreover the horizontal winds having magnitudes smaller than 5m/ 
sec at 100 km reach amplitudes of 45m/sec at 150 km and then decrease again 
to the order of lOm/sec above 200 km altitude. This wind field must be added 
to the wind structure of the diurnal (1,1,1)- mode in Fig. 7 or 8. It may signifi- 
cantly influence the total horizontal wind structure between 100 and 200 km 
especially at low and medium latitudes. 


CONCLUSION 

We did a numerical calculation of the generation and propagation of two 
important symmetric tidal waves within a realistic thermosphere model. These 
fundamental modes are the diurnal wave (1,1,1) which is described by the Hough 
function (1,-1) at lower altitudes and which degenerates into the spherical func- 
tion P 2 2 at thermospheric altitudes, and the fundamental semidiurnal mode 
(2,2,2) which is described by the Hough function (2,2) within the lower non dis- 
sipative atmosphere and which degenerates into the spherical function P 2 2 at 
thermospheric altitudes. Both wave modes are predominantly generated by the 
solar EUV heat input into the thermosphere. In the case of wave mode (1,1,1), 
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P 2, 2^0' mer 'dicnal wind 2^ c o / zona ^ w *nd ^ c q and vertical wind w^ , 
semidiurnal (2, 2, 2)-mode versus altitude generated by the sum of the three 
erators considered. 
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Figure 12. Total pressure, meridional wind and zonal wind versus colatitude 9 of the 
semidiurnal (2, 2, 2)-mode calculated for four different heights. Solid lines: magnitudes; 
dashed lines: phases. 


the contribution of the wave which is generated below 100 km decays rapidly 
with height and becomes insignificant above 180 km. The electric polarization 
field of the geomagnetic Sq current which maps into the upper part of the 
thermosphere generates via ion-neutral drag a significant fraction of the wind 
components of this wave mode at F2 layer heights. The latitude structure of 
this mode constitutes the pressure bulge of the neutral air at thermospheric 
heights the diurnal component of which peaks at 15h L.T. (Jacchia and Slowey, 
1967). The horizontal wind field of the EUV-generated wave above 200 km is 
nearly identical with the wind system derived by Geisler (1967) and by Kohl and 
King (1967). The electric field induced wind component shift the phase of the 
total wind by about 1 hour to earlier times at 300 km altitude. This phase ad- 
vance was expected from a study of the influence of neutral air winds on the 
ionospheric F2 layer (Kohl et al., 1968). 
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In the case of the semidiurnal mode (2,2,2), electric polarization fields of 
E-layer origin as well as the solar semidiurnal gravitational tidal force appear 
to be negligibly small. However, the wave component generated below 100 km 
can penetrate deep into the thermosphere and remains significant even at 400 km 
altitude. Moreover non-linear momentum coupling from the (1,1,1)- mode into 
the (2,2,2)-mode due to the diurnal variation of the ion-neutral collision number 
becomes important at F layer heights. The wave from below propagating as free 
internal wave at thermospheric heights has a phase delay between its density 
and temperature amplitudes of about 3 hours while density and temperature of 
the EUV-induced wave component peak both at about the same time at 13h L.T. 
Due to the interference between the wave from below, the EUV-induced component 
and the drag induced component the temperature amplitude of the total semi- 
diurnal wave peaks at about 15h L.T. whereas the density has its maximum at 
13h L.T. at F layer heights. This phase delay between temperature and density 
of the semidiurnal mode contributes to the difference between the times of 
maximum of the neutral thermospheric density (determined from satellite drag 
data) and the neutral thermospheric temperature (determined from Thomson 
scattering measurements). 
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